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The problem of estimating the thermal corrections to Casimir and Casimir-Polder interactions 
in systems involving conducting plates has attracted considerable attention in the recent literature 
on dispersion forces. Alternative theoretical models, based on distinct low-frequency extrapola- 
tions of the plates reflection coefficient for transverse electric (TE) modes, provide widely different 
predictions for the magnitude of this correction. In this paper we examine the most widely used 
prescriptions for this reflection coefficient from the point of view of their consistency with the Bohr- 
van Leeuwen theorem of classical statistical physics, stating that at thermal equilibrium transverse 
electromagnetic fields decouple from matter in the classical limit. We find that the theorem is satis- 
fied if and only if the TE reflection coefficient vanishes at zero frequency in the classical limit. This 
criterion appears to rule out some of the models that have been considered recently for describing 
the thermal correction to the Casimir pressure with non-magnetic metallic plates. 
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I. INTRODUCTION 

Dispersion forces, frequently called "van der Waals" or 
"molecular" forces, are long-range electromagnetic forces 
arising from quantum and thermal charge and current 
fluctuations existing in microscopic or macroscopic bod- 
ies at thermal equilibrium. In view of their pervasive role, 
from biology to chemistry, from physics to engineering 
[l| , these weak forces have been the subject of intense the- 
oretical and experimental investigations. A distinctive 
feature of dispersion forces is intimately related to their 
long-range character. In fact, while short-range forces, 
like the exchange electromagnetic interaction, are deter- 
mined by the particular microscopic electronic structure 
of atoms and molecules, long-range forces display a uni- 
versal behavior at large distances. This feature is clearly 
reflected in the famous theory of dispersion forces devel- 
oped long ago by Lishitz Q , where the particular features 
of the bodies participating to the interaction, whether 
atoms or macroscopic bodies, can be fully taken into ac- 
count by means of their macroscopic permittivities or 
polarizabilities. 

After over fifty years, Lifshitz theory still constitutes 
the basic theoretical tool universally used by researchers 
in the field, to interpret the results of modern experi- 
ments on dispersion forces. As an example, we mention 
the numerous recent experiments on the Casimir effect 
(for a recent review, see Q), and the beautiful new ex- 
periments on the Casimir-Polder interaction of a Bose- 
Einstein condensate with a dielectric substrate It is 
important to note that the precision of the most recent 
experiments is of a such a level that in order to correctly 
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interpret the data it is now necessary to take into full 
account a number of small corrections, like temperature 
effects, the effect of surface corrugations, patch effects 
etc. (see Ref. Q for details). In particular, the necessity 
of careful electrostatic calibrations in precision measure- 
ments of the Casimir force has been recently emphasized 

In this paper we focus our attention on the influence 
of material properties of the bodies constituting the sys- 
tem. As remarked above, within Lifshitz theory these 
properties are fully described by the macroscopic per- 
mittivities of the bodies participating in the interaction. 
As it is well known, the latter quantities are complex 
functions of the frequency to of the electromagnetic field, 
as an effect of dispersive and absorptive properties dis- 
played by all real materials. Now, in Lifshitz theory the 
free-energy associated with van der Waals forces between 
condensed bodies is expressed by a sum of terms depend- 
ing on the reflection coefficients of the plates, evaluated 
at (imaginary) Matsubara frequencies ^„ = 2'!TnkBT/h, 
where fc^ is Boltzmann constant, T is the temperature, 
and the discrete index n in the sum runs form zero to in- 
finity. As a result, evaluation of Lifshitz formula requires 
knowledge of the reflection coefhcients over a wide range 
of frequencies, extending from zero up to a few times the 
characteristic frequency ujc = c/(2a) of the system, where 
a is the characteristic separation between the bodies. For 
typical separations involved in present experiments, rang- 
ing from a few tens of nanometers up to a few microns, 
the characteristic frequency falls somewhere from the IR 
part of the spectrum, to the near UV. In order to obtain 
an accurate theoretical prediction for the free-energy, the 
common practice today is to rely on detailed optical data 
of the material used in the experiment. The importance 
of using accurate optical data has been recently empha- 
sized by the authors of Ref. @, where it is shown that 
uncertainty in the optical data may easily result in an un- 
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certainty of several percent in the evaluation of Lifshitz 
formula. One point should however be stressed: for the 
purpose of evaluating the n = term of Lifshitz formula, 
it is always necessary to make an extrapolation of the 
reflection coefhcients to zero frequency. This is a very 
delicate point, indeed, because the result strongly de- 
pends on how the extrapolation is done. Understanding 
what is the correct extrapolation is crucial in particular 
for determining the thermal correction to the free-energy 
in systems involving one or more conductors, because 
this correction is strongly affected by the magnitude of 
the n = Matsubara term for transverse electric (TE) 
polarization. This crucial term is determined by the TE 
reflection coefficients of the plates at zero frequency, and 
as of now different recipes have been proposed in the lit- 
erature for this coefficient, resulting in drastically differ- 
ent predictions for the magnitude of the thermal correc- 
tion. The most popular prescriptions used currently to 
model real metals can be grouped in two classes: Drude- 
like models and plasma-like models. The former models 
[3| are characterized by permittivity functions e(w) dis- 
playing at low frequency the same uj^^ singularity of the 
familiar Drude model for the permittivity of a ohmic con- 
ductor, resulting in a vanishing TE reflection coefficient 
at zero frequency. On the contrary, plasma-like mod- 
els are characterized by the uj~^ singularity displayed by 
the plasma model of IR optics, extrapolated to zero fre- 
quency. The latter class of models includes in particular 
the so-called generalized plasma model considered in Ref. 
0. In this model, ohmic dissipation is not accounted 
for, and only dissipation associated with interband tran- 
sitions of core electron is considered. The w"^ singu- 
larity of plasma-like models entails a non-vanishing TE 
reflection coefficient in the limit of zero frequency, differ- 
ently from Drude-like models. The different behaviors of 
the TE reflection coefficient at zero frequency implied by 
the two classes of models have a dramatic impact on the 
Casimir force in the limit of large separations between 
the plates, since Drude-like models lead to a force that is 
about one-half that implied by plasma-like models. We 
remark also that in this limit plasma-like models give al- 
most the same force as the simple ideal-metal model of 
real conductors [1, H, At separations between the 

plates smaller than one micron, where the Casimir force 
can be measured most accurately, the various prescrip- 
tions imply predicted Casimir forces that differ only by 
a few percents, and it is therefore difficult to distinguish 
between them by an experiment. 

The experimental situation is still unclear, because the 
present accuracy of Casimir-force measurements is still 
not sufficient to detect the temperature correction to the 
Casimir force between two metallic bodies, and indeed 
there are at present several ongoing and planned exper- 



iments to measure it 
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A recent accurate experiment 
using a micromechanical oscillator [l^l appears to favor 
the plasma model prescription, but this claim is not yet 
universally accepted by the community [9;]. We remark 
that as of now there is only one experiment that detected 



the influence of temperature on the Casimir-Polder in- 
teraction of a Bose-Einstein condensate with a dielectric 
substrate [3| . Clarifying this problem is important also in 
view of the many experiments on non-newtonian forces 
at the submicron scale, some of which use metallic sur- 
faces at room temperature [l2l |. making it necessary to 
estimate accurately the contribution of dispersion forces 
that must be subtracted from the observed signal. For 
a lucid description of the problems involved in precision 
Casimir experiments we address the reader to the recent 
paper [13]. 

In this paper we shall examine the low frequency pre- 
scriptions for the TE reflection coefficient of insulators 
and conductors, generally used in investigations of dis- 
persion forces, from the point of view of their compati- 
bility with a very well known theorem of classical solid 
state physics, namely the Bohr-van Lceuwcn theorem 
(l^ . This theorem originated early in the 20th century, 
in an attempt to explain the absence of strong diamag- 
netism in normal conductors placed in an external mag- 
netic field. After a simple physical argument due to Bohr, 
who showed that no net diamagnetic currents can arise 
in a bounded conductor subjected to a static magnetic 
field, the theorem was put on a firm theoretical basis by 
H. J. van Leeuwen. In essence the theorem states that 
in classical systems at thermal equilibrium matter de- 
couples from the transverse electromagnetic field. It oc- 
curred to us that perhaps this theorem could be used to 
discriminate between existing models used in the current 
literature to describe dispersion forces. Let us see briefly 
how the connection arises. One observes that the reflec- 
tion coefficients of a surface determine the macroscopic 
response of the surface to an external electromagnetic 
probe placed outside the surface. It is now the essence 
of the famous fluctuation dissipation theorem T5?| that 
such response functions are intimately related to equilib- 
rium averages of suitable macroscopic observables of the 
system. In the case of interest to us, the fluctuation dissi- 
pation theorem relates averages of the fluctuating electro- 
magnetic flelds outside the surface to its reflection coeffi- 
cients [isl ]. This fundamental relation is known to imply 
a set of general constraints, originating from microscopic 
reversibility [l^ , that must be satisfied by the reflection 
coefficients of any real material, like for example impor- 
tant reciprocity relations [13]. Therefore we were led to 
wonder if the Bohr- van Leeuwen theorem can be used to 
put any further constraints on the permitted behavior of 
the reflection coefficients. We remark that the Casimir 
and Casimir-Polder interactions are equilibrium phenom- 
ena, and therefore they must conform to the principles of 
equilibrium statistical physics. In order to put this idea 
to a test we shall evaluate the spectrum of the fiuctuating 
electromagnetic field in the empty space outside one slab, 
and between two plane-parallel slabs, characterized by a 
local dielectric response, carefully separating the longitu- 
dinal and the transverse components of the electromag- 
netic field. Having done this, we shall verify whether 
or not the transverse component of the field decouples 
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from the slab(s) in the classical limit, as required by the 
Bohr-van Leeuwen theorem. Interestingly, we shall see 
that the answer depends exclusively on the behavior of 
the reflection coefhcients of the slab(s), in the limit of 
zero frequency. In this way we obtain a rather stringent 
test to decide whether a definite model is admissible or 
not from the point of view of classical statistical physics. 
The important result is that the Bohr-van Leeuwen theo- 
rem is satisfied if and only if the TE reflection coefficient 
vanishes at zero frequency, in the classical limit. As a re- 
sult, we find that Drude-like models are compatible with 
the Bohr-van Leeuwen theorem, while neither plasma- 
like models nor the ideal-metal model pass the test. Our 
conclusions appear to be consistent with the findings of 
a recent paper which presented a microscopic calcu- 
lation of the Casimir force between two metallic plates, 
in the asymptotic limit of large separations between the 
plates, when the force is dominated by classical thermal 
fluctuations and the theorem is supposed to apply. It was 
found there that in this asymptotic regime, the micro- 
scopic model assumed in [ly] predicts the same Casimir 
force as the Drude prescription, and therefore it is in dis- 
agreement with both plasma-like models and the ideal 
metal results, which we recall both predict in this limit a 
force of double magnitude as that implied by Drude-like 
models. 

We point out that in the recent literature on the ther- 
mal Casimir effect, another criterion based on statistical 
physics has been widely considered to discriminate be- 
tween alternative models for the reflection coefficients of 
a material slab. This other criterion requires that the 
Nernst heat theorem be satisfied, in the limit of zero 
temperature [l^l . This alternative criterion leads to con- 
clusions that are not in agreement with what we found 
on the basis of the Bohr-van Leeuwen criterion, for one 
finds that the Nernst heat theorem is not satisfied by the 
Drude prescription in the idealized case of perfect crystals 
with no defects, but it is satisfied both by the (general- 
ized) plasma models, and by the ideal metal model. In 
the Drude case, Nernst theorem is however restored if an 
arbitrarily small amount of impurities are present in the 
crystal ^] . While we cannot offer a complete resolution 
of this contradiction, we remark that since the Nernst 
theorem is intrinsically a quantum result, this criterion 
is a sense orthogonal to the one proposed in this paper, 
which is essentially classical. In our judgement, in the ab- 
sence of a definitive answer, we note that the Bohr-van 
Leeuwen criterion appears to be more pertinent than the 
Nernst criterion, for a theoretical assessment of the de- 
bated problem of thermal corrections to the Casimir and 
Casimir-Polder effects at room temperatures, since it is 
well known that the difficulties posed by this problem are 
basically of a classical nature (see Sec. V below) . In addi- 
tion, we note that the extrapolation to zero temperature 
of Lifshitz theory poses very non-trivial problems, asso- 
ciated for example with the possible presence of spatial 
non-locality (anomalous skin effect [2^). For recent re- 
views of some theoretical aspects involved by the Nernst 



heat theorem in the context of Casimir physics, we ad- 
dress the reader to Refs. [qI. [23|. 

The paper is organized as follows. In Section II we 
review the fluctuation-dissipation theorem, in the con- 
text of general linear response theory, while in Section 
III the theorem is used to derive expressions for the cor- 
relators of the electromagnetic flelds outside dielectrics 
and conductors. In Section IV we verify if these corre- 
lators satisfy the Bohr- van Leeuwen theorem outside a 
planar slab, for a number of models of dielectrics and 
conductors. The case of a plane-parallel cavity, of the 
type used in Casimir experiments, is considered in Sec- 
tion V. Finally, Section VI contains our conclusions and 
a discussion of the results. 



II. FLUCTUATION-DISSIPATION THEOREM 

In this Section we briefly review the principal results of 
linear response theory, and in particular we present the 
general fluctuation-dissipation theorem for linear dissi- 
pative media. For a review of linear-response theory we 
address the reader to Refs.p^. 

In linear-response theory, one considers a quantum- 
mechanical system, characterized by a (time- 
independent) Hamiltonian Hq, in a state of thermal 
equilibrium described by the density matrix p 

p = e-^" /tr{e~P") , (1) 

where (3 — l/^ksT). The system is then perturbed by 
an external perturbation of the form: 

ffext = - / rf'r ^ Q3 (r, t) f, (r, t) (2) 

j 

where fj{r,t) are external classical forces, and Qj{r,t) 
is the dynamical variable of the system conjugate to the 
force /j(r, t). One may assume without loss of generality 
that, in the absence of external forces, the equilibrium 
values of the quantities Qj (r, t) all vanish: {Qj (r, t)) = 0. 
The presence of the external forces causes a deviation 
S{Qi{r, t)) of the expectation values of Qj(r, t) from their 
equilibrium values. If the forces fj{r,t) are sufficiently 
weak, S{Qi{r, t)) can be taken to be linear functionals of 
the applied forces fj (r, t) according to the formula: 

S{Q,{r,t)) I '^'^ / dt'4>,,{v,v\t~t')f,{v\t') . 

j J J —00 

(3) 

The above Equation assumes that the system was in equi- 
librium at t = —00, and that it reacts to the external 
force in a causal way. The quantities 4>ij{v, r', t — t') are 
called response functions of the system. In principle, they 
can be measured by applying to the system of interest 
suitable external classical probes. We now define the ad- 
mittance (pij (r, r', w) as the (one-sided) Fourier transform 
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of the response function 0.^ (r, r', t): 



(4) 



Being the one-sided transform of a real quantity, the ad- 
mittance satisfies distinctive analyticity and reaUty prop- 
erties. First, it is an analytic function of the complex 
frequency w in the upper complex plane C"*" = {w ^ 
uj + iS, S > 0}. Second, it satisfies in C+ the following 
reality condition 



,-(r,r',-u;*) = C-(r, r', 



(5) 



The latter property in particular implies that the ad- 
mittance is real along the imaginary frequency axis 

By a straightforward computation in time-dependent 
perturbation theory one may prove that the response 
functions 0ij(r, r',i — t') are related to the equilibrium 
expectation values of the commutators of the dynamical 
variables Qi{r,t), in the absence of external forces: 

0,,(r,r',i-i') - A,, (r,r',t-t')^(i-i')- (6) 

Here 9{x) is Heaviside step function {9{x) — 1 for x > 0, 
9{x) = for X < 0) and 



admittance. Consider first the case when Ay (r,r',t) is 
odd in t. This is the case of interest to us, because the 
commutators of the electromagnetic fields are indeed odd 
in time. Then we find: 



dt Ay (r, r', t) e"^* = 2i dt Ay (r, r', t) sin(wt) 



2ilm[0,y (r,r',u;)] 



(10) 



Upon substituting the expression on the second line of 
Eq. (fTO| in the l.h.s. of Eq. ([9]), and after performing the 
inverse Fourier transform of both sides, we then obtain 
the following important relation: 

({Q,(r,t) g,(r',0)}) 



duj 



(11) 



By further exploiting the fact that, by virtue of Eq. 
the imaginary part of the admittance is an odd function 
oi Lu, we can rewrite the above Equation as: 

({Q.(r,t)g,(r',0)}) 



Ay(r,r',t-i')-^([Q.(r,i),Qj(r',t')]> , (7) 



with Qi{r,t) the Heisenberg operator: 



(8) 



As it is well known, starting from Eq. ([6]) it is possi- 
ble to derive several general fluctuation-dissipation theo- 
rems, that allow to express the (symmetrized) correlation 
functions of the quantities Qi{r,t) in terms of the dissi- 
pative component of the response functions (f)ij . The 
form of the fluctuation-dissipation theorem of interest to 
us is expressed by the following relation 



/oc 
diAy (r,r',i) e'- 



I OJ 



dt{{Q,{r,t) Q,(r',0)}) 



(9) 



where 



Ei3(uj) = — coth 



flUJ 

2kBT 



is the average free-energy of a quantum oscillator with 
frequency uj in equilibrium at temperature T, and 
{AB} — [AB + BA)/2 denotes the symmetrized product 
of the operators A and B. If Aij(r, r',t) has a definite 
parity under inversion of time, we can easily verify that 
the l.h.s. of Eq. ^ can be expressed in terms of the 



2 duj 

= - —Ep{uj) Im [0y (r, r', lo)] cos{ujt) . (12) 

Even if we shall not need it, for completeness we report 
also the analogous relation that holds if Aij(r, r',i) is 
even in t: 

({Q,(r,t) g,(r',0)}) 



TT Jo 



duj 



E0{uj)Re[^,j{r,v\uj)]sin{ujt) . (13) 



The above two relations constitute the content of the 
fluctuation-dissipation theorem. Note in particular that 
in the odd case the two-times correlation functions of the 
quantities Qi{r,t) in Eq. p2)) is expressed in terms of 
the dissipative part of the admittance. It is important to 
remark that the integrands on the r.h.s. of Eqs. and 
psp have no singularity ai uj — 0, despite the presence 
of the singular factor oj^^. This so because, from the 
definition Eq. Q, we see that for vanishing frequency 
the admittance 0y (r,r',O) is a finite and real quantity. 
This ensures that the integrands in Eqs. (fT^ and (fT^]) 
both have finite limits as lu approaches zero. 

In view of a next use, it is now useful to derive a 
formula for the equal-time correlators of the quantities 
Qi{r, t), in the classical limit. This can be easily done by 
setting t = and taking the limit ?i ^ in the r.h.s. of 
Eqs. and (fT5|) . The only non-trivial case to consider 
is the odd one, for in the even case we see from Eq. ([Te 
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that the equal-time correlators vanish identically. Then, 
from Eq. we obtain: 

lim({Q,(r,0)Q,(r',0)}) 

h—>-Q 

= ^Im r^0,(r,r',.). (14) 
Jo ^ 

Assuming that the admittance 0ij(r, r', w) vanishes suffi- 
ciently fast at complex infinity in C"*", we can take advan- 
tage of analyticity in C+ of (pij (r , r' , w) to evaluate the 
integral on the r.h.s., by rotating the integration contour 
from the real axis to the imaginary one. The rotated 
contour of integration T consists of an infinitesimal arc 
surrounding the origin in the right sector of C"*" , followed 
by the whole imaginary axis. Since the admittance is 
real along the imaginary axis, it is clear that the part of 
the integral over F extending along the imaginary axis 
does not contribute to the r.h.s. of Eg. p^ . and there- 
fore we find that the imaginary part of the integral results 
entirely from the contribution of the infinitesimal arc sur- 
rounding the origin. After easy evaluation of the latter 
contribution, we obtain the simple result: 

lim({Q,(r,0)Qj(r',0)}) = fc^T^y (r, r', 0) . (15) 

We then reach the important conclusion that in the clas- 
sical limit the equal-time correlators for the quantities 
(5i(r) are simply proportional to the zero-frequency limit 
of the admittance. It is opportune to remark that in 
the above derivation we have implicitly assumed that the 
admittance itself is independent of h. Obviously when 
quantum effects contribute to the admittance, the quan- 
tity ^y(r, r',0) in the r.h.s. of Eq. must be un- 
derstood as the classical limit of the admittance at zero- 
frequency. 



where U{r,t) and A±(r,t) are, respectively, the scalar 
and the transverse vector potentials for the electromag- 
netic field, while p'^"^*-^ (r, t) and j^^^^*' (r, t) denote, respec- 
tively, external classical distributions of charge and cur- 
rent. Note that j^'^'''(r, t) is assumed to be transverse: 

V-j^"'^=0. (17) 

In what follows we shall assume that the bodies consti- 
tuting the system are non-magnetic {fi = 1) dielectrics 
or conductors with sharp boundaries, characterized each 
by a frequency dependent electric permittivity e{uj). It 
is moreover assumed that the bodies are homogeneous, 
in such a way that the permittivities are constant func- 
tions of the position within the volume occupied by each 
body. Under such conditions, it is shown in Ref.[28j that 
we have two independent sets of response functions for 
the scalar and the vector potential: 

U{r,t) = I dt' I (fv'G{Y,Y',t-~t')p'^'''^^\v',t') , 

•J —CO 'J 

(18) 

Ax(r,t) = i /* dt' [ d3r'Gx(r,r',t-t')-jf *^(r',0 , 

(19) 

where G_L(r, r', t — t') has to be understood as a dyadic 
Green function. Recalling that the commutators of the 
electromagnetic fields are odd functions of time, from 
the general fluctuation-dissipation theorem Eq. (|12p we 
obtain the following expressions for the correlators of the 
electromagnetic fields: 

(Mr,<)[/(r',0)}) 
^ / —EB{uj)lm[G{Y,Y',ijj)]cos{ujt), (20) 



III. CORRELATORS OF THE 
ELECTROMAGNETIC FIELD OUTSIDE 
DIELECTRICS AND CONDUCTORS 

In this Section we shall use the methods described in 
the previous Section to derive formulae for the correla- 
tors of the electromagnetic fields outside dielectrics and 
conductors. In the spirit of linear response theory, this is 
done by placing outside the bodies a suitable distribution 
of classical electric and magnetic dipoles, that work as 
external probes for the electromagnetic field [IS]. In or- 
der to separate the longitudinal and the transverse parts 
of the field, we consider an external Hamiltonian of the 
following form [2^ 

i/ext = j d3r[C/(r,t)p(-*)(r,t)-iA^(r,t).jf *^(r,i)] , 

(16) 



({(7(r,t)Ai,(r',0)}) =0, (21) 
({Ax.(r,t)A^,(r',0)}) 

2 dio 

= -/ — Ep{uj)lui[G^^j(Y,Y',uj)]cos{ujt) , (22) 
Ti" Jo ^ 

where G(r,r',a;) and G_L(r,r',a;) denote the one-sided 
Fourier transforms of the Greens functions: 

/•oo 

G{y,y',oj)= dtG{Y,Y',t)e"^\ (23) 
Jo 

POO 

Gi(r,r',c^) = / dtG±{Y,Y\t)e'^' . (24) 
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The Fourier transforms of the Green's functions can be 
obtained by solving the foUowing field Equations implied 
by macroscopic Maxwell Equations: 

V- [e(r,w) VG] = -47r(5(r-r') , (25) 

{A+ e{Y,uj)u^/^)G^{vy,uj) = -47r5i(r-r') , (26) 
where (5j.(r — r') is the transverse delta- function dyad: 

with k — |k|. The above field equations must be supple- 
mented by standard boundary conditions at the bodies 
interfaces, and must be further subjected to the con- 
ditions required for retarded Green's functions. The 
Green's functions G and Gj_ satisfy a number of gen- 
eral properties, that are consequences of microscopic re- 
versibility and of analyticity and reality properties of the 
permittivity e(a;) of any causal material. For a review of 
these important properties we address the reader to Ref. 
[Hi (and Refs. therein). 

For our purposes, it is convenient to split the Green's 
hmctions, outside the bodies, as sums of an empty-space 
contribution plus a correction arising from the material 
bodies: 

G(r, Y',t-t')^ G(°) (r - r', i - i') + (r, v',t^t') , 

(28) 

and 

G^{v.v'.t-t') = Gf (r-r',t-t') + rf "*^(r,r'.t-t') . 

(29) 

Here, G^"-* and G^-* denote the Green's functions in free 

space, while F^™^*^ and F^"''*'' describe the effects re- 
sulting from the presence of the bodies. Such a split- 
ting presents the advantage that all singularities arc in- 
cluded in the free parts G^"^ and G^'' , while the quanti- 
ties F*^"*^*^ and F^^'**'' are smooth ordinary functions of 
r and r'. Upon using Eqs. P5|) and (|^ into the r.h.s. 
of Eqs. (|20ll22p we obtain the following equations for 
the changes of the field correlators outside the material 
bodies, arising from the presence of the bodies: 

5{{U{v,t)U{v\Q)}) 

Jo ^ 

(30) 

<5({t/(r,t)Ai,(r',0)}) =0, (31) 
5{{A^,{r,t)A^,{r',0)}) 

= 1 E^icu) Im )(r, r', u;)] cos{ujt) , (32) 



with an obvious meaning for the symbols. It is now 
easy to derive from the above formulae the expressions 
for the corresponding changes of the equal-time correla- 
tors of the electromagnetic fields that we shall need in 
the following Sections. For the longitudinal electric field 
E|| = -V[/, from Eq. ^ we obtain 

<5({i?p(r,0)i?||,(r',0)}) 

For the transverse electric field and for the magnetic 
field B, since Ex ~ -cr^ dA±/dt, and B = V x Aj_, 
from Eq. ([5^ we obtain: 

5({Fx,(r,0)Fx,(r',0)}) 
^1 r^Epiu;)kllm[F^-'\r,r',co)], (34) 

Jo ^ 

where ko — uj/c and 

Sm{r,0)B,{r',0)}) 

= - / — £;^(c.) Im [ (Vr X X S7r'h] ■ (35) 

TI" Jo ^ 

Finally, from Eq. (fT5|) we obtain the following equations 
for the matter contributions to the equal-time electric 
and magnetic correlators in the classical limit: 

hm 5{{EUr, 0) ^^.(r', 0)}) = keTS^^ljir, r') , (36) 
and 

hm <5({Bx.(r,0)Bi,(r',0)}) = fcs TBit](r, r') , (37) 

where we defined 

£g;(r,r')=lhn(fc^/^:f (r,r',c.)), (38) 

and 

(r, r') = lim ( Vr x x Vr' ).. ■ (39) 

According to the Bohr-van Leeuwen theorem [l3| the 
transverse electromagnetic fields decouples from matter 
in the classical limit. From Eqs. ([36|) and ([37| we see 
that the theorem is fulfilled if and only if the quantities 

^it] ^^AAj ^h vanish identically outside the bodies. 
Equipped with these formulae, we are ready now to ex- 
amine whether the theorem is satisfied in the simple case 
of a dielectric slab. 
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IV. THE BOHR- VAN LEEUWEN THEOREM 
FOR ONE DIELECTRIC OR CONDUCTING 
SLAB 

In this Section we use the results of the previous Sec- 
tions to verify whether the Bohr- van Leeuwen theorem is 
satisfied outside a plane-paraUel dielectric or conducting 
slab characterized by a spatially local permittivity e{Lu). 
We choose our cartesian coordinate system in such a way 
that the z axis is perpendicular to the slab surface, with 
the slab occupying the z < half-space. 

The relevant Green's functions for this problem have 
already been worked out in Ref . [1^ . The matter contri- 
bution f||"^''^ to the tensor Green's functions was found 
to have the following form: 

In this equation, U^J"^^^"^ is the quantity 



We can now evaluate the quantities £^-- and B]^-j defined 
in Eqs. ([38| and ([39|) . Using the following relations 

e^^) ^ Tik^*^ +0(t^') , (47) 

k, =ik^+0{uj^) , (49) 

and observing that k^^'^^ x = ^^'^'> and ^^^^ x k(='=' 
fcg ej_ it is easy to verify that 

+ [Ap\lJ, ki_) - f{Lu)] k^'^k\-^ } g»k<+'-r-*k<-)r' ^ (gg) 



,(P) 



ik<+>T-ik<"'T' 



(41) 



where kj^ denotes the projection of the wave-vector onto 
the (x, y) plane, and we have defined k^ = yjk^ — k'\^ (the 
square root is defined such that lm{kz) > 0), = zxkx, 
k(±) = k_L±/c^z and = fej^z^fc^ki, while r^'\uj ,\<i±) 
andr(p)(w,ki) are the familiar Fresnel reflections coeffi- 
cients for TE and TM waves, respectively: 



e{uj) k^- s 
e(w) kz + s ' 



(42) 



(43) 



where s = t{Lo)kQ — fc^, and again the square root is 

defined such that Im(s) > 0. For the quantity V^^™''' we 
have 

where 



kl dx^dx'j 



27rfc 



ik' + '-r-ik'^'r' 



, z>0, (45) 



where k^^^ = k^ ± iki_ z and the reflection coefficient 
f(ijj) is 



r(w) 



e(w) + 1 



-fcgr(P)(w,fc_L)ej 



-ik^ 



(51) 



As we see from the above two Equations, whether the 
Bohr-van Leeuwen theorem is satisfied or not depends 
entirely on the behavior of the reflection coefficients, or 
what is the same of the electric permittivity e(a;) of the 
slab, in the limit of zero frequency. We consider now 
several classes of models for the low-frequency behavior 
of the permittivity of the materials that are usually used 
in present experiments on dispersion forces 3]. First we 
have insulators, whose permittivities approach a finite 
limit at zero frequency: 



e(w) = £0 + 0{uj) (insulator) 



(52) 



Then we have normal (i.e. non superconducting) non- 
magnetic ohmic conductors. For these materials several 
distinct models have been considered in the current liter- 
ature on dispersion forces, and we classify them generally 
as Drude-like models, plasma-like models and the ideal 
metal model. By Drude-like models we mean any models 
displaying at low frequency the same singular behavior of 
(44) the familiar Drude- model, characterized by an sin- 



gularity: 



0(1); (Drude - like models) , (53) 



where ctq is the dc conductivity. By plasma-like models 
we mean instead any models characterized by the same 
Lo^'^ singularity displayed by the familiar plasma-model 
of IR optics 



(46) e{u;) 



n'i 



■ + 0{uj ^) (plasma — like models) , (54) 
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where fJp is the plasma frequency. These models include 
in particular the so-called generalized plasma model con- 
sidered recently in connection with the Casimir effect [1] . 
Finally we have the ideal- metal model which is better for- 
mulated directly in terms of the reflection coefficients: 

r(P) = f = -A'^ = 1 (ideal metal) . (55) 

We now estimate the quantities S^lj and i3^'^ using the 
above models. It is a simple matter to verify that for 
all models, the reflection coefficients A^'' , r*^*^ and f are 
finite in the limit cu 0. It is also possible to check 
(see also Appendix B of Ref. p^) that the difference 
AP\LU,k±) — f{uj) occurring in the r.h.s. of Eq. ([50|) . 
vanishes always at zero frequency (in the ideal case it is 
identically zero). In view of this, we see from Eq. (|50p 

that the quantities £j_^j vanish for all models: 

£^ = . (56) 

However, in the case of we see from Eq. (IFTI) that 
only the second term between the curly brackets vanishes 
always in the limit of zero frequency, leaving us with: 

=|^,(^)(0,fc^)fcWfc(-)e»fc<+'--'fc'-'r' 

(57) 

where we set 

r(") (0, ki_) = lim A"') (uj, fc^) . (58) 

This equation shows that the key quantity to consider is 
the zero- frequency limit of the TE reflection coefficient, 
for we see that the Bohr-van Leeuwen is satisfied if and 
only if r''*-' (0, A;_l) vanishes. Now, it is easily seen that 
the insulator model and Drude-like models both imply a 
vanishing value for r('*)(0, k±_): 

A^\0, kj_) = (insul. and Drude — like models), 

(59) 

On the contrary, with plasma-like models we find 



For the ideal-metal model we obviously have 

r'-"^ (0, A:_l) = -1 (ideal metal) . (61) 

In view of these formulae, we reach the important con- 
clusion that both the dielectric and the ohmic (or Drude) 
models are consistent with the Bohr-van Leeuwen theo- 
rem, while both plasma-like models and the ideal model 
are not. It is interesting to remark that even when the 

(cl) 

latter two models are considered, the quantity B]_/j be- 
comes negligible at large distances from the slab, and 
therefore the inconsistency revealed here is expected to 
be important only in the study of proximity effects like 
the Casimir effect to be considered in the next Section. 

To avoid misunderstandings, we should repeat the 
warning made at the end Section II. We are tacitly ad- 
mitting here that the reflection coefficient A^^O,k±) is 
independent of h. Obviously, when quantum effects are 
important, it is understood that in Eq. (|57p one should 
take the classical limit of r('^)(0, k±). This consideration 
applies for example both to magnetic materials, and es- 
pecially so to superconductors |24| . 



V. THE BOHR- VAN LEEUWEN THEOREM 
AND THE THERMAL CASIMIR EFFECT 

In this Section, we shall discuss the consequences of the 
Bohr- van Leeuwen theorem for the much debated prob- 
lem of the thermal Casimir effect in metallic systems. 
Therefore, we consider the Casimir apparatus consisting 
of two non-magnetic plane-parallel slabs, separated by a 
vacuum gap of width d. We let eiiuj) and £2(1^) the per- 
mittivities of the two slabs, which are assumed to occupy, 
respectively, the regions z < and z > d. 

The Green's functions for this problem were worked 
out in Ref. [2^, and we can use them here. The mate- 
rial contribution p'^^'*^'' to the transverse Green's function 
was found to be of the form: 

p(cav) ^ -Q(cav) ^ -yicav) ^ /g2) 



r(^) (0, k^) = + "f Z (plasma - hke models) . 

k^ + ^^k\ + Qp/c'^ 

(60) Here 



,ik< + '-(r-r') , zk(-'-(r-r')^ 



'1 

As 



ik'+'-r-ik'^'-r 



As 



ik'"' ■r-ik<+' ■r'+2ik^d 



e±ie±j 



A-2 



. An 



(« 



(+)^(+)g»k<+).(r- 



. ^(-)^(-)gik(-'-(r-r')^ 



'1 
Aq) 



-ik 



(-)r 
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(p) 



(63) 



where 



„(") 



a = s,p are the Fresnel reflection coef- 



feients of slab i for polarization a, and Aa — 1 

1 52vj,(cav) 



rJ^Vj"^ exp(2ifc2(i). As to V''^''^' it has the expression: 



with 

^(cav) _ 



T/(cav) 



27rfc 



(64) 



and F^^^"^ given above, after some elementary algebra 
one obtains the following expressions for (tIP^) and 

;rTi(mat)\ 



duj 



X Im 



„2k I d 



rir2 



(71) 



r2 e 



ik(-'- 



ik' + >-r'-2fc I d 



(65) 

where ^ = 1 — ri{Lo) r2(w) exp(— 2fc^ d). Concerning the 
matter contribution to the scalar Green's function F^'^^'^\ 
in Ref. [2S| it was shown that 



F 



(cav) _ ^(cav) 



(66) 



As it is well known 2] the Casimir pressure acting on the 
slabs is given by the thermal average of the zz component 
^jnjmat)^ of the matter contribution to the Maxwell stress- 
tensor. Since, according to Eq. (PT|) the scalar potential 

and the vector potential are uncorrelated, {Tzf^*"^) de- 
composes into the sum of a longitudinal and a transverse 
contributions: 



(mat) \ 



{T. 



(mat) 



(67) 



and it is interesting to evaluate the two contributions 
separately. Upon recalling the classical expression of the 
Maxwell stress tensor in vacuum: 



1 /I 1 

Tij — ^ I -^OijEkEk — Ei Ej + -SijBkBk — Bi B^ 



we find 



(68) 



(^ifr") ^ j^Y.^k5{{E\\,{v,Q)E\\,{v,Q)]) , (69) 



{Ttf) = r-E,{u:)lu,[r^{u)] (72) 



where 



dfcx k± Iq 



-2ik^d 



-,2k, d 



- 1 



rir2 



(73) 



with q = —ikz- We note that when we take the sum 
of and (T'^™^'''), the longitudinal contribution 

C^ll™*^'^) cancels against the transverse contribution re- 
sulting from the second term between the curly brackets 
on the r.h.s. of Eq. ([73|) . The resulting expression for 
^j^jmat)^ then reproduces the following well known Lif- 
shitz formula for the Casimir pressure P{d, T), expressed 
as an integral over real frequencies 2]: 

P{d,T)^~— —Et,{u;) dk^k^ 
Jo ^ Jo 



X Im 



(74) 



and 



^) = h^^^ ['5({i?i-.(r, 0) E^u{v, 0)}) 



6{{Bk{v,Q)Bk{v,Q)})] , 



(70) 



where Ai 
and 



A2 = -A3 = 1. By recalling Eqs. 
and using the explicit expression of Ft'^^^) 



We can now verify whether the transverse contribution 

(r|™^*^) vanishes in the classical limit, as required by the 
Bohr-van Leeuwen theorem. By taking the limit of the 
r.h.s. of Eq. ([711) for ^ ^ we obtain: 



lim(r('"^')) = -^Im 



Jo ^ 



(75) 



Evaluation of the integral on the r.h.s. is made easy af- 
ter we note that the quantity 71 (w) vanishes for large 
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frequencies, and is an analytic function of the frequency 
in the upper complex half-plane C"*", as a result of an- 
alyticity properties of the reflection coefficients. These 
properties permit us to rotate the contour of integra- 
tion from the real axis towards the imaginary axis, as we 
did already in Section II. Since along the imaginary axis 
71 (w) is real (because the reflection coefficients are real 
for imaginary frequencies) , one finds that the integral on 
the r.h.s. of Eq. ([75)1 receives its only contribution from 
the pole in the origin. The latter is easily evaluated, 
giving us the result: 



— ksT Van TAuj) 



(76) 



We then reach again the conclusion that validity of the 
Bohr- van Leeuwen theorem depends on the behavior of 
the reflection coefficients, or what is the same of the per- 
mittivities, in the limit of zero frequency. At this point 
we assume for simplicity that the slabs are made of the 
same material, in such a way that £1(0;) = e2{ijj) = e('j-'), 
and we consider again the four models Eqs. (|52j|55p for 
the permittivity in the low- frequency limit. It is easy to 
check that for all these models the quantity: 



(r(p)) 



1 



1 



(77) 



occurring on the r.h.s. of Eq. (|73p. vanishes like w^, 
and therefore it does not contribute to 71 in the limit 
of vanishing frequencies. In view of this, only TE-modes 
may give a contribution to 71 in the limit of vanishing 
frequencies and from Eq. fTS]) we obtain: 



2tt 



dk \ k , 



„2k±d 



- 1 



(78) 



The Bohr-van Leeuwen theorem requires that the quan- 
tity on the r.h.s. vanishes for all separations, and this is 
only possible if A^^{0, k±) vanishes. Recalling the values 
for TE reflection coefficient r'-'^^O, fcj_) fisted in Eqs. 
[6T|) . we see again that only the insulator and Drude-like 
models are consistent with the Bohr- van Leeuwen the- 
orem, while plasma-like models and the ideal model are 
not. 

As it is well known after rotation of the frequency do- 
main of integration from the real axis to the imaginary 
axis, Lifshitz formula Eq. ([7^ takes the form of a sum 
over so-called imaginary Matsubara frequencies: 



p(d,r) 



ksT 



E 



1 



dk± k± q„ 



E 



r(")(i^„,fc^)rr(*Cn,fc±) 



- 1 



(79) 



where ^„ = 2TTnkBT /Ti are the Matsubara frequencies 
and qn — Cn/c^- By comparing the r.h.s. of 

Eq. with the r.h.s. of Lifshitz formula Eq. we 
see that the former quantity coincides with the n = 
contribution to the Casimir pressure for TE polariza- 
tion. As it is well known this very term has been the 
object of a long debate in the recent literature on the 
thermal Casimir effect, and as of now there is no consen- 
sus among experts on its actual magnitude in the case 
of normal metallic plates [1, [13] ■ The analysis carried 
here shows that the Bohr-van Leeuwen theorem requires 
that this term be zero for normal non magnetic conduc- 
tors. Quantum effects permit of course non vanishing 
values for the n = Q TE contribution, without violating 
the theorem, for example in magnetic materials and es- 
pecially in superconductors [13] (see remarks at the end 
of the previous Section). 



VI. CONCLUSIONS AND DISCUSSION 

The problem of the thermal correction to the Casimir 
and Casimir-Polder interactions in systems involving nor- 
mal ohmic conductors has attracted considerable atten- 
tion in the recent literature on dispersion forces. De- 
spite numerous theoretical and experimental investiga- 
tions, the resolution of this problem is not clear yet. 

On the theoretical side, several distinct models have 
been proposed, that give widely different predictions for 
the magnitude of this correction. The difficulty is that 
in order to get a definite value for the thermal correc- 
tion, it is necessary to make a definite extrapolation to 
zero frequency of the optical data of the plates, because 
the theoretical prediction is very sensitive to the limit- 
ing behavior of the reflection coefficients for vanishing 
frequency. Unfortunately, this extrapolation cannot be 
done solely on the basis of direct measurements, and it 
unavoidably requires making some theoretical assump- 
tions. Different ansatz have been proposed in the litera- 
ture, each supported by some arguments, that give widely 
different predictions for the magnitude of the thermal 
correction. 

The experimental situation is not definite either. In 
principle, the best way to clarify the issue would be to 
measure the Casimir force between two metallic plates at 
separations of a few microns, because for such large sep- 
arations different theoretical models predict forces that 
differ by as much as about fifty percent. Unfortunately, 
at these distances the Casimir force is very small, and 
at the moment no one has been able to measure it with 
sufficient precision [ll|. Below one micron, where the 
Casimir force can be measured most accurately, the ther- 
mal correction gets very small and therefore it is difficult 
to observe it. At the moment of this writing, the most 
relevant experiments in this regard are those at Purdue 
University [l^ (see also Ref.[^). While the achieved pre- 
cision is not sufficient to detect the thermal correction 
predicted by certain theoretical models (of the plasma 
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type), the authors claimed that their measurements rule 
out at high confidence level the much larger thermal cor- 
rection predicted by alternative theoretical models (of 
ohmic conductor type). As of now, no other experimen- 
tal groups have carried out measurements of the thermal 
correction to the Casimir force. For a definitive assess- 
ment, it would be desirable to have more experiments, 
possibly of different types. Perhaps, further insights may 
come from superconductors [23 |. 

In this paper we have examined the most widely used 
models for the low-frequency behavior of the reflection 
coefficients of dielectrics and conductors, from the point 
of view of their consistency with the Bohr- van Leeuwen 
theorem [l3|. This is a well known theorem in classical 
statistical physics, stating that at thermal equilibrium 
the transverse electromagnetic field decouples from mat- 
ter. As it is well known, this theorem provides the ba- 
sic explanation for the absence of strong diamagnetism 
in normal conductors. The theorem has a very gen- 
eral character, and it fact it also holds under the less 
restrictive assumption of local kinetic equilibrium [25| . 
In this paper, we evaluated the correlation functions for 
the transverse electromagnetic field outside a dielectric 
or conducting slab, and inside a planar cavity of the type 
used in Casimir experiments. Upon taking the classical 
limit of these correlators, we found that decoupling of 
the transverse electromagnetic field occurs if and only 
if the reflection coefficient for transverse electric (TE) 
modes vanishes at zero frequency, in the classical limit. 
According to this result, we conclude that the dielectric 
and the ohmic conductor models are completely consis- 
tent with the Bohr- van Leeuwen theorem, while neither 
plasma-like models nor the ideal metal model are. In- 
terestingly, in the case of a cavity, the average value of 
the Maxwell stress tensor, providing the Casimir pres- 
sure, is consistent with the theorem if and only if the 
n = Matsubara mode for TE polarization gives no con- 
tribution to the Casimir pressure. This term is precisely 
the source of the present controversies on the thermal 
Casimir effect for metallic plates. From the point of view 
of the Bohr- van Leeuwen theorem the answer is clear: 
this term must be zero, up to quantum effects. We re- 
mark that non-vanishing values for the TE reflection co- 
efficient are of course possible in materials displaying a 
magnetic response, and in superconductors, since both 
phenomena arise from quantum effects that disappear in 
the classical limit. When such materials are considered 



a on- vanishing contribution to the Casimir pressure from 
the n — Matsubara mode is clearly possible, without 
violating the Bohr- van Leeuwen theorem. 

The reader may be disturbed by hearing that the fa- 
miliar ideal-metal model is inconsistent with the Bohr- 
van Leeuwen theorem. After all, this is the model that 
Casimir himself used to investigate the effect that goes 
under his name, sixty years ago [26j . We remark that 
this model constitutes an extreme idealization of real 
metals, and therefore one should not expect it to be 
universally valid. A hint that this model may provide 
incorrect answers in problems involving proximity ef- 
fects of the quantized electromagnetic field near material 
surfaces was pointed out long ago by Milonni [S^ who 
found that ideal-metal boundary conditions imply a vio- 
lation of equal-time canonical commutation relations for 
the transverse electromagnetic field in the vicinity of a 
metallic surface. In a recent work [28] , the author of the 
present paper showed that such a violation is an artifact 
of the idealized boundary conditions for a perfect metal. 
Indeed, if account is taken of dispersive and dissipative 
features of real materials, no such violation is found, and 
the canonical commutation relations are restored at all 
point outside the conductor. As to the range of valid- 
ity of the ideal-metal model, we have seen above that 
the field correlators derived from this model violate the 
Bohr- van Leeuwen theorem only at finite temperature 
and only in the vicinity of the metallic surface. We infer 
from this that the model has to be considered unreliable 
in the study of the proximity effects like the Casimir and 
the Casimir-Polder interactions, in situations where tem- 
perature effects become important. 

Finally, we note that in the recent literature on the 
Casimir effect there has been a lively debate on the in- 
fluence of free charge carriers that always exist in poor 
conductors and in semiconductors [1, . The Bohr- van 
Leeuwen theorem cannot give any help to discriminate 
between the alternative models that have been proposed 
for these materials, because in this case differences in the 
predicted magnitudes of the Casimir and Casimir-Polder 
interactions arise either from longitudinal fields or from 
TM fields. As far as longitudinal fields are concerned, 
the Bohr- van Leeuwen theorem obviously does not ap- 
ply. In the case of TM electromagnetic fields, we have 
seen in Sections IV and V that they always decouple in 
the classical limit, whatever model is used for the dielec- 
tric function. 
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